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Abstract 

Arcade-type magnetic field structures originating from tfie photosphere are relevant to the un- 
derstanding of different types of solar prominences, coronal loops and coronal heating. In most of 
the existing literature, these loop-like magnetic structures are modeled as force-free fields (FFF) 
without any plasma flow. The system is assumed to be isolated and non-dissipative. In reality 
the photospheric plasma hardly qualifies to be isolated in nature and the existence of an external 
drive in the form of a dynamo field is always possible. Solar prominences are not ideal either since 
dissipative effects are believed to play a major role in coronal heating. The above observations 
indicate that a force-free model based on a non-dissipative plasma may not be a suitable candidate 
to replicate the arcade structures and further investigations are required. In our work, we have 
proposed the arcade structures as minimum dissipative relaxed states (including both the viscous 
and resistive channels) pertaining to a two-fluid description of the plasma. The obtained relaxed 
state is non force-free in nature and appropriate to an open system with external drives. The 
Euler-Lagrange equations are solved in Cartesian coordinates subject to the existing photospheric 
boundary conditions. The solutions are seen to support flow-containing arcade like magnetic field 
configurations with inherent dissipative properties that may play an important role in the coronal 
heating. An interesting feature observed is the generation of different types of arcades with the 
variation of a single parameter characterizing the relaxed state. Also, recent observations with the 
LASCO coronagraph on board the SOHO spacecraft suggest that the helmet streamers originating 
from the sun may have an internal triple-arcade structure. The two-fluid relaxed state obtained 
here is also seen to support such structures. 
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I. INTRODUCTION 



d line(s) joining two magnetically active 
^. The active regions themselves are 



Arcade type magnetic structures play an essential role in solar physics, particularly re- 
garding the coronal loops, coronal heating and coronal mass ejection. A typical coronal loop 
consists of a bunch or a single arch-like magnetic fie 
regions of opposite polarity on the photosphere 
separated by a neutral line, where the normal component of the magnetic field is zero. In 
general the coronal loops are observed to be long-lived. By long-lived we mean the life-time 
of the structure involved is more than that expected from a conventional magnetohydro dy- 
namic (MHD) analysis. The coronal loops are seen to have life-times varying from an hour 
to more than one day whereas the MHD time scale for a typical loop parameter is of the 
order of seconds j4| only. 

The global structure of the magnetic arcades are modelled either by assuming a potential 
(current-free) field Q] or by considering the plasma to be force-free jsl Q, US ^^^^ the 
magnetic field satisfying the equation 

V X B = AB (1) 

with A as a spatial function or a constant. The fff model is justified by an an order of 
magnitude analysis that shows the structures to be of low /3_|l0| or on arguments similar 
to that of Taylor in the context of laboratory plasmas Il2i|. where it is envisioned as 
a relaxed state characterized by the magnetic energy minimum along with an invariant 
magnetic helicity 0, Q] • The second scheme appears to be more lucrative, since by definition 
it makes the arcade structures long-lived and at the same time imparts a sense of stability to 
them. This stems from the fact that the relaxed state being a spontaneously achieved state, 
any deviation from it is expected to generate the same mechanisms that were responsible to 
attain the relaxed state at the very first place. As a result, the system will once again try 
to converge to its earlier state and the perturbations will die out. In other words, there is 
always an anticipation that the relaxed states are also inherently stable by nature. But this 
is only an expectation and a detail analysis is always required to study the relevant 

stability properties. 

Although the force-free description is successful in explaining the gross topological fea- 
tures of the magnetic arcades but at the same time also suffers from major drawbacks. As 
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observed through EUV spectrum P| the plasma temperature has a monotonically increasing 
profile inside the loop be.ng minimum at the axis and maxtaum at the surface. The density 

shows no significant variation transverse to the loop axis jj] and even can be considered to 
be uniform 15]. These observations imply that the pressure is not an uniform parameter 
inside the loop and a force-free description may not be an appropriate candidate to model 
the arcade. In a recent work by G. A. Gary it has been recognized that the solar corona 
itself has regions of different plasma (3. It is only the mid-corona sandwiched between the 
photosphere at the base of the coronal loop and the upper-corona where the plasma (3 << 1 
and the force-free condition can be assumed to hold good. The other two regions are essen- 
tially characterized by > 1 and hence support non force-free fields. Based on the SXT 
Limb Data 17], Gary has proposed a model that admits a range of j3 values in the mid 
corona. Also in an alternative scenario the solar fiares are pro pos ed to be natural conse- 
quence of the high /3 plasma confined in a curved coronal loop 18|- The above-mentioned 
works essentially points that a force-free model is not adequate in explaining all the features 
of arcade fields and a non force-free description is highly desirable. 

Another feature missing in the force-free description of the arcade fields is the plasma fiow. 
In all the force-free models mentioned above, the plasma is assumed to be ideal or nearly 
ideal. This essentially decouples the magnetic and fiow fields, making them independent of 
each other. In contrast to this finite plasma fiows have been observed in the coronal loops 
jlfll ]. a typical form being the Evershed fiow which is directed out of the sunspots at the 
photosphere and into them in the chromosphere Q]. A finite pressure difference between 
the foot points is believed to be responsible for such flows. Flows are important in context 
of the solar winds also where the plasma is seen to flow outwards along the open fleld lines. 
The importance of flow has recently been recognized by Mahajan et. al. j2ll] also. According 
to them viscous dissipation does play an important role as long as the flow fleld is treated 
as an integral part of the plasma dynamics and accounts for the coronal heating during the 
very formation of the loops. 

From the above discussions it emerges that although the force-free model is suitable in 
explaining the gross topological features of the arcade-type magnetic flelds, further studies 
must be made in order to accommodate features such as non-uniform pressure proflle and 
flnite plasma flow. To include such features, in the following we propose the arcade flelds as 
driven relaxed states governed by the principle of minimum dissipation rate (MDR). Before 
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going any further, we want to declare that our goal is rather modest and is confined only 
to the feasibility study of the arcade structures as MDR relaxed state. In this paper we 
will be more concerned with the general topological features and leave all the quantitative 
analysis; like the heating rate or the magnitude of the plasma flow for future works. Also 
for analytical progress, essential in understanding the natural coupling of the self-consistent 
magnetic flelds and non-zero plasma flows, we will assume the density to be constant. The 
situations where such an assumption is valid being already been discussed by Mahajan et. 



al 



2l|, no further clariflcations are given here. 



II. FORMULATION OF THE PROBLEM 

To formulate a variational problem for obtaining the relaxed state, we need a minimizer 
(minimization of which yields the relaxed state) and invariants or constraints that essen- 
tially determines how far the system will be propped away from its trivial outcome. A 
crucial point to be noted is, most of the arcade flelds, including the solar corona, are exter- 
nally driven magnetoplasma systems. For the solar corona, the external drive is provided 
by the convective motions in the bounding photospheric surface penetrated by the coronal 



magnetic fleld 



22j |. With dissipative processes playing a key role |21|], such systems are not 
in thermodynamic equilibrium and we have to look for suitable minimizers and invariants. 
Principle of minimum dissipation rate flrst proposed by Montgomery and Philips in the con- 



text of plasma physics 



23| is one of the suitable candidates to describe such off-equilibrium 



m\M (and 



magnetoplasma systems related to various fusion oriented conflgurations 
the references therein). We must point out that in absence of a rigorous mathematical proof, 
principle of MDR still enjoys the status of a conjecture- the validity of which depends on 
the agreement of predictions based on it with the experimental observations. The success 
of MDR in predicting experimentally realizable relaxed states in laboratory plasmas makes 
it a viable candidate for the astrophysical plasmas also. 

In the following we describe the plasma by two-fluid equations. This has two basic 
advantages over its single-fluid or MHD counterpart. In two-fluid description, at-least in 
context of the relaxation physics, plasma flow gets coupled to the magnetic fleld in a natural 
way. Also, the two-fluid description being more general than MHD, it is expected to capture 
certain results that are otherwise unattainable from a single-fluid description. In two-fluid 
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formalism the magnetic helicity is replaced by the generalized helicities defined as 



K^= Pa- n^dr (2) 
Jv 

where Pq, and are the canonical momentum and canonical vorticity for the a species and 
are given by 



I 9o » 

c 

(3) 

ri„ = V X p„ 

As can easily be acknowledged, the expression for in equation (2) is gauge independent 
only if the boundary condition fla ■ n = is satisfied. Such a boundary condition makes the 
system to be isolated from the surroundings. To accommodate arcade structures, which are 
inherently driven systems j^], we must modify the above definition of generalized helicity. 
An alternative expression for the generalized helicity can be found by extending the definition 



to the two-fiuid case. Let us define 



of Jensen and Chu 



Ka = J^Pa- n^dT -J^K- Kdr (4) 

where Pq, and P^ are different inside the volume of interest and are same outside. The 
gauge invariance can easily be seen by writing Pq, = Pq + Vtp, ip being some arbitrary 
scalar continuous at the boundary. 

K'^ = K^ + ftPin^-n'J -Ma (5) 

where the surface integral is to be evaluated over the bounding surface S of the volume V. 
In writing the above equation we have used the Gauss' theorem along with the fact that 
V ■ Q,a{^a) = 0- So the generalized helicity Ka as defined by equation (4) is gauge invariant 
provided the boundary condition 

(riQ -n'j-h = o (6) 
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is obeyed. If we further assume that there is no flow-field coupling at the surface i.e., the 
flow and the magnetic fields are two independent variables and the variation of one does not 
affect the other, we have 



(B - B') • n = 

(7) 



{uJa -Uj'J-h=0 



at the surface, where 



UJa = V XUa (8) 

is the vorticity. Equations (7) then represent a type of natural boundary conditions inherent 
to the problem. Next, with generalized helicity given by equation (4), in the following we 

derive the helicity balance equation. 

For the a component the momentum balance equation is given by. 



maUa-^ = riaQa + X B^ - V/la + /^aV^Ua - UaqaV^ (9) 

where we have assumed the flow to be incompressible and ha = Pa + 9, where g is the 
gravitational potential and p is the pressure. Writing E = — V0 — {l/c){dA/dt) and using 
the vector identity 



(Ua ■ V)u« = - U« X (V X u„) (10) 



the momentum balance equation can be cast as. 



9Po 

dt 



= -V 



Ua 2 



+ Ua X fio 



(Vxfu.-M£(Vx)B 
Un 47r 



(11) 



where we have used the notation (Vx)" = V x V x ...to n terms. Taking curl on both sides, 
we get the vorticity balance equation as 



_^^ = VxKxf2„)-^(Vx)^u„ 
at Tin 



qavc 

An 



(Vx)2B 



(12) 
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Using the above two equations along with the definition (4) for the generahzed hehcity, 
corresponding balance equation for the a component is obtained as, 



a 




dt 



(13) 



with 



niaU' 



2 1 




Ga — 



An 




2 



'a 



(14) 



47r Ha 



= Prandtl number 



qacria V 



While deriving the above equation, we have used the boundary condition described by 
equation (6). In the above balance equation the surface integrals represent the helicity 
injection terms while the volume integral represents the dissipation rate for K^. For an 
isolated system characterized hy fl'^, ■ n — 0, Ua ■ n — and n x P = (the last condition is 
due to the fact that in isolated systems with conducting material boundary, n x A — and 
77, X Uj = which is the no-slip boundary condition and is applicable to viscous systems) , only 
the dissipation term in the balance equation survives and for a static plasma characterized 
by Uq = is similar to the dissipation rate of the magnetic helicity. 

To proceed further, in the following we consider a two-component hydrogen plasma. 
Neglecting the electron mass over the ion mass with = —qg — q and the quasineutrality 
to hold good, we have 




qrjc 
An 



(15) 
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J _ 47r /ii 

Li — 

qcn 7] 



In writing the last expression we have neglected the electron viscosity over the ion viscosity 
due to its smallness by the electron to ion mass ratio. 

Following Montgomery et al j^l, let us assume that K a. bounded function. Then the 
time average of K over a reasonably long time will make left side of the above equation a 



constant. Prescribing the surface terms 



29| then it is possible to form a variational problem 



where the relaxed state is obtained by minimizing the total dissipation rates, both ohmic 
and viscous, with the helicity dissipation rates as constraints. In other words, in our model 
we have assumed the field variables to be fixed at the surface. This assumption is more 
of a mathematical requirement to keep the subsequent analysis simple and at per with the 
variational principle in classical mechanics. The minimization integral is then obtained as. 



J 



(V X B) ■ (V X B) + ^^^'^^' ui ■ Ui 



dr 



+ Ai y"[V X (B + LiLUi)] ■ flidr - Xe J x B) ■ n^dr 



(16) 



the first two terms in the above integral represents the ohmic and viscous dissipation rates 
while the second two integrals represent the ion and electron generalized helicity dissipation 
rates. Aj and A2 are the corresponding Lagrange undetermined multipliers. The Euler- 
Lagrange equations are obtained by equating the first order variation of the above integral 
to zero and treating (5B and 5uj as independent variations and are given by, 



(Vx)2b + i^i±^(Vx)B + ^ + ) (Vx)^, = 
c 2 V c 



(Vxfuji + ^^(Vx)lj, 



1 fqLi 



+ mi) {VxfB = 



rriiXiC 2miLi \ c 

Eliminating B (uj) in favor of Uj (B), the EL equations can be rewritten as 



(17) 
(18) 



(Vx)2b + 



(Aj + Xe)q 27mq Aj / qL 



+ 



rajXiC AraiLi \ c 



(Vx)B 



27mq^ Xi + Xe 
H ^ ; B = Vip 



niiC^ Xi 



(19) 
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{Xj + Xe)q ^ 2'Knq 



X, 



mjXjC AmjLj \ c 



qLi 



(Vx)tUi 



27Tnq^ Xi + Xe ^ 
+ ^ ; = Vx 



A,: 



(20) 



where ip and x satisfy Laplace's equation 



V'^ix) = (21) 

A few comments are necessary for the above set of equations. For a constant gauge, the 
equations reduce to the double-curl Beltrami equations obtained by Mahajan and Yoshida 
(8r| . Also the equations are in some sense robust to a variation in the Prandtl number. For 
both the limits Li ^ and Lj — oo and for finite Lj the equations maintain their double- 
curl nature. As the finite Lj case also includes the ideal limit [rj —>■ 0, fii ^ 0), this shows 
that the flow does play an important role in securing a relaxed as well as a steady-state in 
the ideal limit, as conjectured by Montgomery and Philips 



III. SOLUTION OF THE EL EQUATIONS: ARCADE STRUCTURES 

To obtain solutions of the EL equations pertaining to arcade structures let us employ the 
cartesian coordinates. We assume that the x — y plane represents the photospheric surface 
and arcade structures are extended to the positive half z— plane. So our volume of interest 
is the region characterized by the positive values of z coordinate. 

To solve the EL equations, let us define 



2nnq'^ Aj + A^ 
m,c^ A,' 



Bt = B - ; "" v^ 

(22) 



I 27mq'^ Xi + Xe^ 

u)\ =Ui Vx 

rriiC Aj 

In terms of these newly defined variables, the EL equations can be written as. 



(Vx)2Bt + 
2-Knq^ Aj + A 



{Xj + Xe)q ^ 2'nnq 



X, 



rjiiXjC 



AlTLiLi 



qLi 



+ m, 



(Vx)Bt 



A, 



■Rt = 



(23) 
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2 . (Aj + \e)q , 27mq f qLi ^ ^ 



(Vx)W + 



c miXiC 4:miLi \ c 

2 



:vx)LUi 



^2nng_K + K^^ = (24) 

The above equations suggest that the daggered fields can be expressed as a superposition of 
two Chandrasekhar-Kendall (CK) eigenfunctions, i.e., 



(25) 

a;l = Yi + /?Y2 



where Y^'s satisfy the relation, 

V X Y,- = A,Y,- (26) 

The a and P are constants and quantifies the non force-free part in B''' and cuj respectively. 
We assume that the x — y plane represents the photospheric surface and the arcades are 
along the z-direction. In cartesian coordinates a convenient representation for Y with y 
symmetry is 

Y = V X {(piy) + X V X {(piy) (27) 
A 

For (26) to be true, should satisfy the Helmholtz equation 

+ X^(j) = (28) 

With y-symmetry is given by. 



cos{xVW+X^) exp[-kz] (29) 



Since the arcade solutions are bounded only at 2; = and dies out as z ^ 00 jSl, we have 
assumed an exponentially decaying solution along the z-direction. Also, we have written 
only the cos part so that the different field components agree with that given in the works 
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of Browning et. al. [tI]. From equations (27) and (29) different component of the daggered 
fields may be written as, 



and 



where 



Bl 



00. 



UJ 



- \\ cos(/«ia;)e~^v^''i"^? + oi^k\ - cos(K2a;)e"'^v^ 



2 \2 



Aicos(Kix)e V''? + (x\2 cos(K2x)e ^v^l^ 
fti sin(fi;ix)e~^V'^i"^i + afi;2 sin(K2a;)e~^v^^^ 



^/s:f-A?cos(Kia;)e-"v/^ + /^^Z/ti - Ai cos(/€2x)e-"v^ 



2 \2 
2~^2 



Al cos(Kix)e"V''?~^? + /5A2 cos(K2a;)e"V«^i-^i 



Kisin(/tix)e ^v^^ + /3k2 sin(/t2x)e ^v^i^ 



(30) 



(31) 



- + 



J = 1,2 (32) 

The gradient part in equation (22) can be calculated by solving (21) in the region described 
by 2 > 0. 



ip{x) = E{F) sin(/a;)e 



-Iz 



(33) 



Note that we have once gain assumed an exponentially decaying solution along the z-direction 
to be in conformity with the arcade structures. From equation (22), different field compo- 
nents are obtained as, 



By 

B. 



I - Af cos(Kia;)e"V''?-^? + a^/ti - A^ cos(K2x)e~ V'^i-^i + E cos{lx)e~ 
Al cos(Kix)e~^v^''i~^? + aA2 cos(/t2x)e"V'"i-^i 



-Iz 



(34) 



fi;isin(fi;ix)e ^V''? ^? + aK2sm{n2x)e ^V"! ^1 + Esin(/a;)e" 



■Iz 



and 
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\Ik\ - Af cos(Kia;)e"V'^?-^? + /Jy^/ti - Ai cos(K2x)e~V'^i-^i + F cos{lx)e 
Ki sm(Kia;)e"^V''i^^i + (3k2 sm{K2x)e~^y^^^^ + Fsm{lx)e~' 



(35) 



where parameters kj, Xj (for j = 1, 2), a and /3 are to be determined from the boundary 
conditions. In the following, we attempt a simplified solution without any crossings and 
knots. The condition for that is 0, 3ll | 



ki = k2 = k = I 



or equivalently 



(36) 



"1 ~ "-2 ^^2 

The field components for this case then reduce to 



(37) 



cos(kix) + a cos(k2x) + 



E 



^4 - A? 



cos X\/ K 



By = [Ai cos(/tix) + aA2 cos(/t2x)] e 



— kz 



5, 



Ki sin(fi;ix) + aK2 sin(K2x) + E sm{xJ k\ — Af 



"1 



-kz 



A? 



and the vorticity components are obtained as, 



-kz 



(38) 



^4 - A? 



cos(Kia:;) + /5 cos(k2x) + 



V'^? - A? 



cos(a;y^ h\ — Af 



-kz 



ijjy = [Ai cos(kix) + PX2 cos(/t2a;)] e" 



-kz 



Ki sin(Kia;) + [3k2 sin(K2a;) + F sin(x\/ /tf — Af 



(39) 



A point to be mentioned is, as the boundary conditions (6) or (7) has to be obeyed at all 
the points on the x — y plane, we have no other option but to assume that this condition is 
satisfied automatically. In other words we are assuming the photospheric fields are to be of 
the same form as the above set of equations with 2; = 0. In the above set for B if we put a = 
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E = 0, the correspon ding magnetic field represents the photospheric boundary conditions 
used by Finn et. al. 32|- Here, we are assuming a more general but structurally similar 
photospheric boundary condition. Equations (6) or (7) then do not yield any condition 
through which different parameters can be determined. To establish the possibility that the 
MDR relaxed state does support meaningful flow-containing non force-free arcade structures, 
in the following we consider the simplest two-dimensional non force-free state described by 
the condition, 



Al — — A2 — A 

For such a choice, different field components can be written as. 



(40) 



where 



E 



= \/7^ - 1 (1 + tt) cos(7x) + / 2 _i cos(xy^72 - 1) 



By={l~ a) cos(72;)e 



(41) 



(1 + 0)7 sin(7x) + E sm{xJ 7^ — 1 



and E, n and A are redefined as 



K 



(42) 



K = kL 

X = XL 



(43) 



L being the characteristic length scale representing the width of the arcade in the x — y plane 
and X and z are normalized w.r.t it. The coefficient a is obtained from the requirement that 
at a; = 1 the magnetic field lines are open. For that we need 



B.,. = B,i 



(44) 



which gives 
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a = 1 



(45) 

^_ 2y/i^ - 1 COS 7 

COS ■\/7^ — 1 

With these values for a and E, different field components are obtained as, 



B,, = 



r 



2 cos 7X 



2 cos 7 



cos VT 



: cos(_xa/7 



(46) 



B. = - 



27 sin 7X — 



2^7^ — 1 cos 7 . 



sm[x\ry 



sin a/7^ — 1 

As we have pointed out earlier, the resulting magnetic field is two-dimensional with By = 0. 
This is just a simplified case and more elaborate three-dimensional magnetic fields can be 
obtained for other choices of the eigenvalues. But, here our objective is rather modest; 
to model and look for arcade-type magnetic fields as possible MDR states. To do that a 
two-dimensional magnetic field would suffice. The third dimension essentially adds twist to 
the magnetic field lines and is believed to be responsible for S'-shaped solar prominences 
0, 0|. We leave this as our future work. The Arcade structures are obtained by solving 
the magnetic field line equation, 



(47) 



dz _ B, 
dx Bj. 

and are depicted in figures (l)-(5) for different 7 values. It is interesting to observe that 
along with the single-arcade structure, the MDR relaxed state also supports multiple-arcade 
solutions. The triple-arcade solutions are particularly interesting as recent observations with 
the LASCO coronagraph on board the SOHO spacecraft suggest that the helmet streamers 



originating from the sun may have an internal structure similar to that 



33 



34I . The observed 



triple-arcade structures existed for several days and occasionally went unstable leading to 
a new and extraordinarily huge kind of coronal mass ejection (CME). Although a direct 
method for predicting CMEs is not included in our model, but following Bellari |3ll| a 
rough qualitative scenario for the same can be visualized. According to Bellan |3l!] the 
loss of equilibrium occurs when the altitudinal decay factor 7 in equations (46) satisfies the 
condition 
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7 < 1 (48) 

For 7 < 1 the magnetic fields no longer decay with increasing height from the solar surface 
but are oscillatory in nature. In the footsteps of Bellan, we may also propose the following. 
As 7 decreases the arcades passes from more stable to unstable states and as the instability 
threshold is crossed, loss of equilibrium occurs with erupting prominences. In other words, 
this model requires the multiple-arcade structures to be more stable compared to their single- 
arcade counterparts and indeed, in a recent work based on resistive MHD, Wiegelmann et 
al have shown that triple streamers are usually more stable than a single streamer j^^. 
Alongside, we want to put a word of caution. The argument presented here is only qualitative 
and indicates that the MDR relaxed state has the potential to incorporate the phenomenon 
of CME. An actual stability calculation for the states described by equations (46) can decide 
the issue. 

Although, equation (20) supports field-aligned flows as a non-unique solution obtained 
with a special gauge, in the following we look for more general solutions. Different flow 
components can be calculated by realizing. 



dz 

dux duz 

dz 
dm. 



dx 



(49) 



dx 



where we have dropped the subscript i for convenience. From the first and the last set of 
the above equations, y component of the flow-field is obtained as. 



[1 + /?) C0S7X + 



cos(xa/72 — 1 



(50) 



Other two flow components are obtained from the second of the equations (49), which gives 



dur du. 



(1-/5) cos7xe-"vV^ 



(51) 



dz dx 

Assuming an exponentially decaying form for the flows, the flow components in x and z 
direction can easily be calculated from the above equation and is obtained as. 
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Wl'^ - 1 - 7) COS 



COS ^^7^ — 1 cos(7x) — C0S7 cos(xy'7^"-^) 



cos -i/72 — 1 sin(7x) — cos 7 sin(xA/ 7^ — 1 



(52) 



COS - 1(^7^ - 1 - 7) 
To obtain (3 and F we need additional boundary conditions on flow, which are not available 
at this time. Actually by demanding Uy to be zero at the edge, one of the constants may be 
calculated but the other remains still unspecified. Also, we have to look about the matter 
that in our calculations we have obtained By = with a non-zero Uy. In solar physics it 
is hypothesized that due to the shear flow (the y-component) the footprints of the arcades 
move on the x — y plane. If so, we must generate a 3D magnetic field with non-zero y- 
component. Then the choice Ai = — A2 will not do and we have to look for more general 
relations, in other words more general boundary conditions are required. 

Figure (6) shows the x-variation of for single-arcade solution with 7 = 1.5. The 
direction flips sign at the summit of the arcade and hence resembles an evershed flow. The 
corresponding flow hues are depicted in figure (7). 

IV. CONCLUSIONS 

In summary, in this work we have proposed a formulation suitable to obtain externally 
driven relaxed states in two-fluid plasma. The final relaxed state is obtained by utilizing the 
principle of minimum dissipation rate (MDR), which is appropriate for such systems. One 
novel feature of this formalism is the presence of non-trivial flow field coupling that is absent 
in the corresponding single-fluid MDR states. The Euler-Lagrange equations obtained are 
double-curl in nature and represent non force-free configuration. For a constant gauge this 
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reduces to the double curl Beltrami equation obtained by Mahajan and Yoshida 
Since the double-curl Beltrami equation represents one of the steady states of the ideal 
plasma, this highlights the possible role of the non-zero plasma flow in securing a steady, as 
well as, MDR relaxed state. 

To establish the plausibility of MDR principle in predicting arcade structures, an attempt 
has been made to solve the Euler-Lagrange equations in terms of Chadrasekhar- Kendall 
eigenfunctions subject to the arcade geometry in cartesian coordinates with axisymmetry. 
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A solution is obtained for the two-dimensional case characterized by By = 0. One interesting 
feature supported by this new relaxed state is the prediction of continuous transition from 
the single arcades to multiple arcade type solutions with increase in the eigenvalue 7. A 
particular case of interest is the observation of triple-arcade structures as the MDR relaxed 
state. The importance of such structures can be realized from the recent LASCO observa- 



tions on SOHO spacecraft that suggest 
have an internal triple-arcade structure 



;he helmet streamers originating from the sun may 
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, 3^. Also a hypothetical picture of CME can be 
postulated where the arcades transit from multiplets to singlets before the instability thresh- 
old 7 = 1 is achieved. This is consistent with the recent findings that the triple streamers are 
more stable than the single ones js^. In addition, the relaxed state predicts self-consistent 
plasma flow, which for the single-arcade solutions resembles the evershed type. 

The above findings definitely points out that the MDR relaxed states applied to astro- 
physical plasmas is a worthy case for further investigations. Any such investigation should 
involve all the three components of the magnetic field. Only then it will be possible to study 
the effects of twist on the field lines. In principle, this should generate knots and further 
studies can be made on the predictability of S'-shaped solar prominences. 
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FIG. 1: Magnetic field lines corresponding to single arcade structures. 7 = 1.5 
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FIG. 2: Magnetic field lines corresponding to intermediate structures between single and double 
arcades. 7 = 3.0 
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FIG. 4: Magnetic field lines corresponding to intermediate structures between single and double 
arcade. 7 = 3.247 
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FIG. 5: Magnetic field lines corresponding to triple arcade structures. 7 = 3.36 
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FIG. 6: Uz profile for 7 = 1.5 
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FIG. 7: Flow lines lines corresponding to single arcade structures. 7 = 1.5 
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